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By Warren A. Tucker 


SUMMARY 


The problem of Increasing the range of supersonic aircraft by the 
use of twisted and cambered wings is considered, pr imar ily for the pur- 
pose of developing a rational method for the selection of a design lift 
coefficient. Relations and curves are presented from which a suitable 
selection may be made, depending on the relative importance of maximum 
range and top speed. 


INTRODUCTION 


The inherently low values of lift-drag ratio and the resultant 
short ranges characteristic of conventional aircraft configurations 
operated at supersonic speeds have stimulated research on more efficient 
shapes for use at these speeds. This research has of necessity been 
directed along two main lines, one the reduction of the drag at zero 
lift and the other the reduction of the drag due to lift. Along the 
latter line, one of the most promising developments has been the use 
of twist and camber to reduce the wing drag due to lift for a given 
plan form. Contributions to this field have been made by Robert T. 

Jones (refs. 1 to 3), E. W. Graham, and his coworkers (refs, li- and 5), 
and S. H. Tsien (ref. 6), among others. 

One of the principal problems with which the reports just mentioned 
are concerned may be stated thus: Given a wing plan form operating at a 

given lift and a given Mach number, find the shape (as expressed by the 
angle of attack, the twist, and the camber) which^will result in the 
lowest drag. Attention Is directed to the words "at a given lift and a 
given Mach number" since the optimum shape will change both with the 
lift and with the Mach number. For any particular aircraft, however, 
practical considerations will usually dictate the use of a wing having 
fixed twist and camber, even though the aircraft may be required to 
operate over a range of lift coefficient as well as a range of Mach 
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number. Out of this situation arises the problem of how best to select 
the one particular shape that will represent the best compromise for all 
flight conditions , which amounts to selecting the lift coefficient at 
which the drag is to be minimized at a given Mach number. The present 
paper offers some suggestions concerning this selection process. 


SYMBOIS 


CD 

C D ( 

C D, 


drag coefficient, 


Drag 


'0,f 


0,w 


C - D m,w 

AC D 

cl 

°IdL 

% 

Vw 

^■^optjf 

^■Copt,w 

^^opt,w 


Dynamic pressure X Reference area 
drag coefficient of flat wing at zero lift 

drag coefficient of particular minimum-drag wing at' 
zero lift 

drag coefficient of particular minimum-drag wing at 
design lift coefficient 

lowest drag coefficient for a particular minimum-drag 
wing 

drag coefficient measured from the flat-wing zero-lift 
drag coefficient, Cj) - „ (used with the various 

subscripts) ’ 


lift coefficient, 


Lift 


Dynamic pressure x Reference area 
design lift coefficient “ 

lift coefficient at which drag of particular minimum-drag 
wing equals drag of flat wing 

lift coefficient at which v is obtained 

lift coefficient at which (L/P) T11flv ^ is obtained 
lift coefficient at which (V®)max,w obtained 


lift coefficient at which 


(l i/2 /d) 

y / / /max,w 


is obtained 
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^ L / D )max,f 

maximum lift-drag ratio for flat wing 

(l l/2 A)aax,f 

maximum (lift)^/^-drag ratio for flat wing 

(k/D)max,w 

maximum lift-drag ratio for particular minimum-drag wing 

( l 1 / 2 / b )mk,v 

maximum (liftJ-^/^-drag ratio for particular ml nimum-H-rag 
wing 

Kf 

drag-rise factor for flat wing, 

K v 

drag- rise factor for family of minimum- drag wings , 

SCd/^Cl 2 

p = yii 2 - l 


M 

Mach number 

m 

cotangent of sweepback angle of leading edge 

= Sc L /^ 



a angle of attack 

Superscripts : 

* condition of maximum 

** condition of maxi muni 
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ANALYSIS AND DISCUSSION 
Factors Influencing the Range 


It is anticipated that aircraft operating at supersonic speeds will 
visually he powered hy a reaction propulsion system (turbojet, ram jet, 
or rocket, for example), for which the fuel rate (weight of fuel used 
per unit time) will he closely proportional to the thrust, rather than 
to the horsepower, as in the case of ' a propeller-driven airplane. For 
such aircraft, the maximum range will he attained hy flying at the con- 
dition of (L/D) TOAy if the flight Mach number is specified, or at the 


condition of 


( Ll / 2 / D )max 


if the ratio of the atmospheric pressure to 


the wing loading is specified (specifying this ratio is roughly equivalent 
to specifying the altitude). These relations are illustrated in figure 1. 


The characteristics of the propulsion system enter into the pre- 
ceding statements only to the extent that the fuel rate is assumed to 
he proportional to the thrust (for example,, the changes in available 
thrust with changes In altitude or Mach nvuriber do not affect the con- 
siderations) . In the case of any specific aircraft, however, the engine 
characteristics will he a primary factor in the determination of the 
actual operating point for maximum range. Although a detailed discus- 
sion of this problem is not the purpose of this paper, it is mentioned 
so as to emphasize the fact that the operating point for a particular 

aircraft may he neither at (L/D) TT1SY nor at Q?-f^/'P) mar} hut at some 

point between these two conditions. Therefore, after the main assump- 
tions and relations have been stated, the following analysis is divided 
into two main parts, one dealing with (L/P) 771a y and the other dealing 


with 
covered. 


( Ll/£ / D )max' 


so that the limits of probable design conditions are 


Basic Assumptions and Relations 

The drag coefficient of the flat wing is assumed to he expressible 
as a parabolic function of the lift coefficient for a given Mach number 
and wing plan form: 


CD - C Do ^ f + KfC L 2 


( 1 ) 
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This equation is plotted in figure 2. As long as the relation between 
Cj) and Cl is parabolic, the curve can be regarded as applying either 
to a complete aircraft configuration or to an isolated wing. For con- 
venience, however, the term flat wing" is used herein when referring 
to the curve . The drag of the family of wings which at ary value of Cl 
has the lowest possible drag is assumed to be known (from the work of 
ref. 6, for example, or from some other suitable source) and to be given 
by the following equation: 


C D = Cu 0>f + K/Jl 2 (2) 

This equation is also plotted in figure 2. 

The first step is to write the equation for the drag of one of the 
family of minimum-drag wings, that is, of a wing which has a fixed twist 
and camber. Now the addition of small, amounts of twist and camber to a 
flat wing introduces increments of drag and lift which do not vary with 
angle of attack. The drag of a wing with fixed twist and camber can 
therefore be represented by the following equation: 


C D " C Dm,v + *f( °L - c Im,w f 


( 3 ) 


where C-p^ w and Cj^ ^ are shown in figure 2. If equation (3) is to 

express the drag of one of the family of minimum-drag wings, the value 
of Cd given by equation (3) must first be equated to the value given 
by equation (2) for some particular lift coefficient C]^, which is 

called the design lift coefficient. In addition, since the drag given 
by equation (2) is the lowest possible at any value of Cl, then the 
curves represented by equations (2) and (3) must be made tangent at 
Cl = Cj^. These two conditions serve to determine Cj^ v and Ci^.^ 

in equation (3), and the equation for the drag, of the particular 
minimum-drag wing becomes 


- %, f + - l )*- 0 ^ 2 + K * 


Cl- 
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00 


The relations Just discussed are illustrated in figure 2. The discus- 
sion of the curve represented by equation 00, in particular its rela- 
tion to equation (l), forms the remainder of the paper. Some miscel- 
laneous relations for the various coefficients are presented in 
appendix A. 
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Design Considerations Based on (L/D ^y 

Maximum lift -drag ratio .- As mentioned previously, twist and camber 
may be used to produce greater values of the lift-drag ratio, and in 
particular of the maximum lift -drag ratio (L/D) TT1fty w , than can be 
obtained from the flat wing. From equations (4) and (l) the ratio of the 
maximum lift-drag ratio for any particular minimum-drag wing (i*/^)max,w 

to the corresponding value for the flat wing (L/D) max p can be obtained 
as a function of the two quantities K^/Kp and f> wiiere 

Cx^pt p is the lift coefficient corresponding to (L/D) mflv p. The ratio 
is as follows: 


(^A^max.w 

( L / I) )max,f 



For the flat wing the following relations hold: 


( 5 ) 


(L/D) 


max,f 



C 


■^opt ,f 



> 


( 6 ) 


The ratio given by equation ( 5 ) is plotted in figure 5(a) • Of some 
practical interest is the fact that for a given value of K^/Kp (which 

amounts to a given plan form and' Mach number) there is a particular design 
lift coefficient at which the greatest increase in (L/D) max is obtained. 

(l/d) 

• * max w 

If an asterisk is used to denote the condition of maximum *— , 

( L / D )max,f 

then the following relation can be written: 


f(L/D) 


T* 


max,v 


(L/D) 


max,f 


% ,1 
CL opt,f 


( 7 ) 
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Equation ( 7 ) is shown as the dashed line in figure 3(a)* One criterion 
is thus provided for the selection of C]^, if the only item of concern 

in the particular design under consideration is the max im um possible 
increase in (L/D)^, to the exclusion of other items (such as increases 

in minimum drag) . 

The numerical values of the (L/D)^ ratio in figure 3(a) are also 
of interest. As shown in appendix B , values of K^/Kf in the neighbor- 
hood of 0.5 or 0.6 can probably be expected for reasonable plan forms and 
Mach numbers. From figure 3(a), increases in (L/D) maY of 30 to 40 per- 
cent can therefore probably be realized by the introduction of the proper 
twist and camber. 


Other characteristics .- In connection with the lift-drag ratio, the 
lift coefficient corresponding to (L/D) ffiax is of importance. The ratio 

of CLqpt v to ^Lopt f ts given by the following equation: 



This ratio is, of course, always greater than unity. 


In the case of configurations for which the maximum speed capabili- 
ties, as well as the maximum range, are important, the increase in mini- 
mum drag caused by the addition of twist and camber to the flat wing must 
be considered. (Although a detailed study should include the effect of 
Mach number, the change in minimum drag with Mach number is small enough 
so that for the purposes of this paper the constant Mach number case can 
be considered.) This increase in minimum drag can be expressed as 
follows : 


^z - 1 + Vi M / Cl * V 
° D o,f " J 


(9) 


and the lift coefficient at which the min imum drag occurs can be deter- 
mined from 


= L _ KyN ° L d 

CL opt,f V L opt,f 


( 10 ) 
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If, 

maximum 


as before, 
< L / D >max,f' 


an asterisk is used to denote the condition of 
then equations (8) to (10) become 


c w \* = CL d = 1 

C Lopt,f 1 CL opt,f 


(11a) 



(lib) 


(The reason for the seemingly strange result given by this equation 
for the impractical (or ideal) case of Ky. — » 0 can be understood if a 

plot is made similar to fig. 3(c) but with Ci^/Cx^pt w aB the abscissa. 

The case of Kf > » can be visualized easily by making a sketch similar 

to fig. 2.) 


Equations (8) to (ll) are plotted in figures 3(b) to 3(d) . It is instruc- 
tive to compare the maximum increase in (l/D) ^ r v with the corresponding 

increase in minimum drag, that is, the values given by the dashed lines 
of figures 3(a) and 3(c), which come from equations (7) and (lib). These 
values are platted in figure 4 and show that for a realistic range of 
K w /Kf (greater than 0.4) the percentage increase in minimum drag is 

always greater than that in (L/D) ^ ^., if the design condition is chosen 
to give the greatest possible increase in (L/D) max . This situation is 

very undesirable, of course, if the top speed of the configuration under 
study is important, since the top speed is greatly dependent on the mini- 
mum drag. 

Evidently if the top speed of the configuration is important, a more 
satisfactory design condition would be one which would yield something 
less than the maximum possible (L/D) TTia y increase but which would keep 


\ c % 1 

C Lopt,f j ° L opt,f 0 Ld/ CL opt,f 'jl\ l /K e 
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the minimum-drag increase within acceptable limits. Such a design con- 
dition is possible, as shown by the data presented in figure 5- As an 
example, consider a curve with K^/K ^ = 0 . 5 . The maximum possible 

(L/D) max increase of 4l percent is achieved only at the expense of a 

50-percent increase in minimum drag. However, if the design condition 
is modified so that a (say) 30-percent increase in (L/D) TnpY is obtained 

(by choosing a value of c L d / c L opt; f of about 0.68), then the minimum- 

drag penalty is decreased to 12 percent. Similar considerations apply, 
of course, to any of the other curves with K^/Kf constant. As indi- 
cated in the example just discussed, the lines of constant C^^Lop-^ p 

give the proper value of design lift coefficient to use once a satis- 
factory design condition, as measured by the gain in (L/D) ^^ balanced 

against the penalty in minimum drag, is selected. 


Design Considerations Based on 




Considerations analogous to those developed in the preceding section 
entitled "Design Considerations Based on (l/p)^^" can be expressed on 

the basis of (L 1 / 2 /p) Triay . This development proceeds in a manner similar 
to that for the case of ( L/d) mpY . 


Maximum value of L ' /D .- The equation for 
sponding to equation (5) is 

( l 1 /£ >Lc,w 

k 1/2 /°W 


( l1/2 / d )« 


corre- 


4 - aef - (■ - 8 2 er - 4 - 


For the flat wing, the following relations hold: 


, 1 / 2 ' 
J L 

C D , 


max,f ^C^// 2 


1 

’’•^opt , f 


2m 

% 
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It is emphasized that C]^^. ^ has the same meaning here that it had in 
the previous development; that is, it is the lift coefficient corresponding 

to (L/D) ( not ( L^/d) A 

' 'max,f \ ' ' 'max,f/ 

The ratio given by equation (12) is plotted in figure 6(a). If a 

(l l/g / P)max,w 

(l 1/2 /b) ' 


double asterisk is used to denote the condition of maximum 
then the following relation holds: 


(l 1//2 /d) 

' ' 'max,w 

- 

% . 

1 


(l/d) 

max,w 

(l 1//2 /d) 

' / max,f 

CL opt,f 



(l/d) 

max,f 


max,f 


(14) 


Equation (14) is plotted as the dashed line in figure 6(a) and provides 
a criterion for the selection of Cj^, the basis of selection in this 

case being the development of the maximum possible increase in 
, to the exclusion of other factors. 


max 


Other characteristics .- The lift coefficient corresponding to 

\L ' /d) 0 „ „ Is designated herein by Ct ’ . The ratio of Ct* , 

' max,w * -Mopt,w n>pt,w 

to Ct is given by the following equation, which corresponds to 


opt , f 

equation (8) for the case of (L/D) 


max’ 


CI w i % 

— 

1 ^ + J ? +j 

\ _ M 

AM 

C Lopt , f ^ C;i ^opt , f 

_ * l(%/°Bopt,f ) 2 1 

Kfj 

\ 


(15) 


The equations for minimum drag and the lift coefficient at which minimum 
drag occurs (eqs. (9) and (10)) apply to the present case as well as to 
the previous one . 

If, as in equation (l4), a double asterisk is used to indicate the 


condition of maximum , . - 

(l 1 / 2 / d) 

take the following forms: 


— x - >- , then equations (15), ( 9 ), and (10) 


max,f 
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opt ,w 
-'opt , f 


GL d _ 1 = 1 / opt,v 

C;L opt , f 'fiXw/Kf ) f i 


/ C - D m,v ' 

1 “ 4 1 


' k 

v%tJ 

3 K'Wt.f) 

2 3< 

\ Kfj 


c % 1 ill _ = w^ 

'opt fij ° L opt,f 3CL d/ CL opt,f l^Kw/Kf 


Equations ( 9 ), (10), (15)* and- (l6 ) are plotted in figures 6(b) to 6(d). 

A comparison of the maximum increase in U 1/2 A>Lx ■with the corre- 
sponding increase in Cp^ is given in figure 7- The values in figure 7 

come from equations (l^) and (l6b). The comparison is rather more favor- 
able than the corresponding comparison in figure b-, but in some cases the 

balance between ( l1/2 A>Lcv and Cjj^ v represented by figure 7 may 

be inacceptable . In such cases a more suitable design compromise can be 
selected from the curves of figure 8, which corresponds to figure 5 Tor 
the (L/D) max case. Considerations similar to those discussed in con- 
nection with figure 5 apply To figure 8. 

Comparison of Designs Based on ( ■ L/P ) TnpY and on ( L"*"^/D) max 

Because of the existence of at least two design bases (that is, 

the question naturally arises as to what 

extent a design which emphasizes one of these factors is penalized with 
respect to the other. Some indication may be had from the answer to the 
following specific question and its converse: If ^ is chosen 

to give ^ L//D ^maX;y | } -what percentage of will be 

UA>W.f 
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realized? The answers to these questions are shown in figure 9> together 
with the minimum- drag increase and the design lift coefficient for each 
design condition. The first point to he made from the figure is that, 


a i/2 A u 


is 


regardless of which factor ^that is, (L/D) mav or 

maximized, the increase in the other is always at least 92 percent of the 
maxi mu m possible increase. The second point is that the minimum-drag 


increase resulting from maximizing is much less than that 

resulting from maximizing (L/D) max . Therefore, a design lift coeffi- 
cient chosen to give the maximum possible increase in ( L^/^/P ) max will 

result in a wing which will develop at least 93 percent of the maximum 
possible increase in (L/P) mRV but which has a much smaller minimum drag 

than if the design had been chosen to give the maximum possible increase 
in (L/P) mRY . Such a design lift coefficient therefore suggests itself 

as a good choice as long as the minimum drag is not too critical. For 
cases in which closer attention must be paid to the minimum drag, a more 
suitable compromise can be selected from figures 5 and 8. 


CONCLUPING REMARKS 


Under the assumption that, for a given wing plan form, the lowest 
drag obtainable by the use of twist and camber is known, information has 
been presented which will allow the rational selection of a design lift 
coefficient. It is suggested that, for a given Mach number and wing plan 
form, a design lift coefficient which gives the largest possible increase 

in the maximum (lift)"''/ -drag ratio is a reasonable choice as long as the 
minimum drag is not too critical. If the increase in minimum drag due to 
twist and camber must be kept very small, a somewhat lower design lift 
coefficient should be used. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September 8, 195^* 
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APPENDIX A 

MISCELLANEOUS RELATIONS FOR THE VARIOUS COEFFICIENTS 


Certain relations, although not directly hearing on the main purpose 
of the paper, are nevertheless of incidental interest and are easily 
derived. For example, if the values of Qd given by equations (l) 

and (4) are equated and the resultant lift coefficient is called Ct 
( see fig. 2), the following simple result is obtained: 


C Le _ 1 


(Al) 


Other similar relations can be obtained. For example, if ACq is 

used to denote the drag coefficient measured from the flat -wing zero -lift 
drag coefficient as a reference, that is. 


= C D - C Do ^ f (A2) 

then these relations are as follows (see fig. 2 for the physical meaning 
of the various coefficients): 


^•m,v _ ^Pm.,w ■, 
c L d AC Dd Kf 


^^m,w _ Ky 
AC DOjV = Kf 


ACn 


0,w 


ag d , 



(A3) 


The preceding relations are shown in figure 10. 
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APPENDIX B 


PROBABLE RANGE OF Ky/K f 

In reference 2, Robert T. Jones gave an" equation for the minimum 
drag (at a fixed lift) of a slender wing lying near the center of the 
Mach cone. This equation was also given later by Adams and Sears in 
reference 7> with some discussion of the derivation. For wings with 
zero tip chord, the equation may be written in the notation of the pres- 
ent paper as follows : 


1 v _ 1 1 + 2(ftm) 2 

P^ A/m itgm 


(Bl) 


The drag for the corresponding flat wing, if the leading edge is assumed 
to develop no thrust, is simply 


so that 

*w 
% A/m 

Values of can be obtained from various sources , such as 

reference 8. 

Values of K w /Kf calculated from equation (B3) are plotted In fig- 
ure 11. Also included are some more exact values for the conical trian- 
gular wing (A/m = 4-), taken from reference 6. The agreement between 
these values and those calculated from equation (B3) is good. Although 
equation (Bl) is not strictly applicable to a wing with A/m = 6, the 
values in figure 11 are included to show the trend of the variation of 
Ky/Kf with A/m. The lower values of Ky/Kf obtained for the higher 

values of A/m indicate that the use of proper twist and camber should 
be more beneficial for arrow-type plan forms than for diamond-type plan 
forms. Values of Ky/Kf in the neighborhood of 0.5 or 0.6 can probably 
be expected for reasonable plan forms and Mach numbers. 

As mentioned previously, the values of Ky/Kf presented in fig- 
ure 11 are based on the assumption that no leading-edge thrust is devel- 
oped in the flat .wing. For Increasing amounts of leading-edge thrust, 
the values will move closer to unity. 


i Kf = 1 

P 1 




PCn l + 2(pm)‘ 


(B2) 


/-n T \ 
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* ■* 4 

• • I i I i ii 1 Ml i I ! ' " 1 i i ' ! | I 


MCA TN 3317 



MCA TN 3317 


r 



F/at wmg 
Cp~ * KfC L 




n 

'"l-Q ^IT}} 14 '' 



Q 


Figure 2.- Basic relations. 
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Figure 3 
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Figure 6 
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(a) Design lif t coefficient chosen to give (b) Design lift coefficient chosen to give 

maximum possible increase in (L/D)^. maximum possible increase in (l 1 / 2 /^^. 


/ . x 

Figure 9.- Comparison of increases in (L/D)^^ and and corresponding increases in Cj^, 

for two design conditions. 
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